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Abstract 

The HQET/NRQCD Lagrangian is computed to order ajrnJ. The computa¬ 
tion is performed using dimensional regularization to regulate the ultraviolet 
and infrared divergences. The results are consistent with reparametrization 
invariance to order 1/m^. Some subtleties in the matching conditions for 
NRQCD are discussed. 
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I. INTRODUCTION 


Heavy quark effective theory (HQET) |l[] and non-relativistic QCD (NRQCD) are 
two effective theories that describe the interactions of almost on-shell heavy quarks. HQET 
describes the interactions of quarks of mass m in which the momentum transfer p is much 
smaller than m. The HQET Lagrangian has an expansion in powers of p/m. HQET is 
typically applied to hadrons containing a single heavy quark, such as the B meson, in which 
p ~ Aqcd, the scale of the strong interactions. The HQET expansion is thus an expansion 
in powers of Aqcd/^^- NRQCD describes the interactions of non-relativistic quarks, and is 
typically applied to QQ bound states such as the T. The NRQCD Lagrangian also has an 
expansion in powers of 1/m. The momentum transfer in NRQCD is of order mn, so that the 
small expansion parameter in NRQCD is the velocity v. The size of a term in the NRQCD 
Lagrangian can be estimated using velocity counting rules [Q . The basic difference between 
HQET and NRQCD can be seen from the hrst two terms in the effective Lagrangian, 

c = qUiD^)q + Q^q. ( 1 ) 

^ ^ 2m 

In HQET, the hrst term is of order Aqcd, and the second term is of order Aq^j-j/m, whereas 
in NRQCD both terms are of order mn^. As a resnlt, the quark propagator in HQET is 
-1- ie), and in NRQCD it is 

_^_. ( 2 ) 

{hP — k2/2m -|- ie) 

The HQET/NRQCD Lagrangian is computed in this paper to one loop and order 1/m^. 
Only the terms bilinear in fermions are considered here. There are also fonr-qnark operators 
in the effective Lagrangian. Their coefficients are order and can be obtained simply from 
tree-level matching. 

II. MATCHING CONDITIONS AND POWER COUNTING 

The HQET effective theory matching compntation is a straightforward generalization 
of known results to order 1/m^ One can compnte diagrams in the fnll and effective 
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theories, and match to a given order in 1/m. Since the HQET propagator is m independent, 
the HQET power connting is manifest — one connts powers of 1/m directly from the vertex 
factors. This means that graphs with a vertex of order 1/m'’ do not make any contribntions 
to terms of order 1/m®, with s < r to any order in the loop expansion. 

The nse of NRQCD as an effective field theory is more snbtle. NRQCD with the propaga¬ 
tor Eq. (^) cannot be nsed as an effective Lagrangian to compnte matching corrections, since 
the velocity power connting breaks down.Q The matching conditions for NRQCD should be 
computed using the HQET power counting, by expanding in p^/m. After the HQET La¬ 
grangian has been computed, it can be used for computing bound state properties using 
the NRQCD velocity power counting rules. In other words, the NRQCD propagator Eq. 
should be thought of as the inhnite series 


1 k2 

rTT +-- 


1 


( 3 ) 


2m 


(P — k2/2m -|- ie) 


where one uses the right hand side inside any ultraviolet divergent Feynman graph. This is 
necessary when a cutoff such as dimensional regularization is used to regulate the Feynman 
graphs. NRQED matching conditions have previously been computed using a momentum 
space cutoff . In this case, there is no difference between using the left or right hand sides 
of Eq. (Q). However, a momentum space cutoff can not be used in NRQCD, since it breaks 
gauge invariance. 

The difference between using the two forms of Eq. in a loop graph can be illustrated 
by a simple example. Consider the integral 



A typical NRQCD loop integral has the form Eq. (^). The power a increases as one considers 
more and more divergent loop graphs in the effective theory. The denominator of a typical 


would like to thank M. Luke for extensive discussions on this point. Some of these issues will 


be discussed in a futnre publication. See also |^|. 
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NRQCD loop graph has poles at ~ where p is the external momentum, and at ~ m? 
where m is the quark mass. Thus the scales mi and m 2 in Eq. can be taken to be of 
order p and m, respectively. One can see immediately that the NRQCD power counting 
breaks down. Loop graphs with insertions of higher dimension operators are divergent, and 
can be proportional to positive powers of m because of the term in Eq. (p. The 

positive powers of m from the loop integral can compensate for inverse powers of m in the 
coefficient, and the entire effective Lagrangian expansion breaks down. Now consider the 
same integral, but hrst expand 


1 


i m^ ^ 


+ m2 k'^ k^ 
evaluate the integral, and then resum the series. The answer is 

TT {miY 


dk^ 


(5) 


+ ml) [k'^ + ml) sin tto mf — ml 

The integral is missing the term since it is non-analytic at the origin for a 7 ^ integer, 

which is where the integral is evaluated in 4 — e dimensions. Equation (^) only has inverse 
powers of the high momentum scale m 2 ~ m, and leads to an acceptable effective held 
theory. Thus NRQCD and HQET matching conditions are computed in the same way, and 
the two Lagrangians are the same. 


III. THE LAGRANGIAN 


The continuum NRQED effective Lagrangian at one-loop has previously been com¬ 
puted 0 using a photon mass to regulated the infrared divergences, and a momentum 
space cutoff. This procedure cannot be used in a non-abelian gauge theory such as QCD. 
The kinetic terms in the NRQCD Lagrangian at one-loop have previously been computed 
by Morningstar IjTOI using a lattice regulator. The computations in this letter will be done in 
the continuum using dimensional regularization for the infrared and ultraviolet divergences, 
and for on-shell external states. This has the advantage that one can freely use the equa- 
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tions of motion to reduce the number of operators in the effective Lagrangian . The most 
general effective Lagrangian to order 1/m? (up to held redehnitions) is 


^ f n cr • B 

C = QHiD'^ + C2 — + ca + Cpg —— 

/m im 


cd9- 


[D-E] 


+icsg 


TCp'p g 


cr ■ (D X E - E X Dl 


Sm? 

cr-DBD + DBcrD 


cwi g 




— cw2g- 


D* cr B D* 


8m^ Am? 

D [D X Bl + [D X Bl D 


8m3 


icm g- 


+cai g 


B2 - £2 
8m^ 


ca 2 g 


8m^ 

2 E2 . orr . / E2 \ 


+icBi g" 


<T ■ (B X B - E X E) 

8m^ 


iCB2 g 


cr ■ (E X E) 

8m^ 


>Q, 


( 6 ) 


(7) 


which is the HQET/NRQCD Lagrangian in the special frame v = (1,0, 0,0), and the no¬ 
tation of has been used. The covariant derivative is + igA^°‘T°‘ = {D^, —D). 

Covariant derivatives in square brackets act only on the helds within the brackets. The 
other covariant derivatives act on all helds to the right. The subscripts F, S and D stand 
for Fermi, spin-orbit, and Darwin, respectively. The last seven terms in Eq. (P) are not given 
in Ref. 0, since they were not required for the computation done there. The last four terms 
can be omitted for QED. 

In an arbitrary frame, Eq. (P) can be written as 


Qv ^ ^71 ' V C 2 ~ “t” C 4 “ r 

2m 8m'^ 


cpg 


o'agG 

Am 


ap 


Cog 


D/G 


afS 


vxCT^0{Dl,G^f^} {Dl,a^pG^^} 

+^cs g -^ + cwi g- — 777 -^^—^ - cw 2 g 


8m^ 
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8m^ lOm^ 
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8m^ 


( 8 ) 


D/ GnfiD 


2 ^ol(3^ I ^ ^2 ^ ^ol 


16m^ 


+ ca 2 g 


8m^ 

G„,„G>^?VaVi3 


±'^ag-^^±X j . D±a 

— icM g - 




D 




16m^ 


-F ca 3 g^ Tr 


-icbi g 






16m^ 


*cb 2 g 


<^0/3 


G^a, 




lOm^ 


m r^aB' 
16m^ 

Qv) 


+ Ca 4 5'^ Tr 


8m^ 

' GpaG^^^Vc,Vp 


16m^ 


where 


D': = -v'^ V ■ D. 


(9) 
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The tree level matching conditions can be obtained by integrating out the antiquark com¬ 
ponents, and making a held redehnition to eliminate terms with v ■ D acting on the quark 
helds. The “standard” form of the HQET Lagrangian after integrating out the antiquark 
helds is 


C 


V 


Qy\iv-D + ilp±- - Qv 

L 2m + tv ■ D ) 

Qv ^iv- D - -^Ip^Ip±_ + ■ D)Ip^ 


8m^ 


-?*± (“«' D 


( 10 ) 


The held redehnition 


_ D1 {iv ■ D) D±a gv\D±aG°‘^ 

8m^ 16m^ 16m^ Ibm^ 

. a^pDl {ivD)Di _ . gvy,ao,pDlG^\ ^ 

16 m^ Ibm^ -I ^ 


( 11 ) 


(where the a matrices are understood to be PycrPy) can be used to eliminate the time 
derivative terms, and put the Lagrangian into the “NRQCD” form. The result is Eq. (H) 
with C 2 = Ci = cf = cd = cs = cwi = cai = cbi = 1, and Cw 2 = Cp/p = cm = ca 2 = = 

ca 4 = cb 2 = 0. The ca and cb terms are quadratic in the held strengths, and are order g'^. 
The one loop corrections to these terms will not be computed here. 


IV. QUARK FORM FACTORS AND MATCHING CONDITIONS 


A loop diagram in QCD is a function F{{p}, m, /r, e) where {p} are the external momenta, 
m is the quark mass, /i is the scale parameter of dimensional regularization,and the computa¬ 
tion is done in d = 4 —e dimensions. As an example, consider the diagram Fig. 2, which gives 
radiative corrections to the form factors Fi (g^) and F 2 (g^). In dimensional regularization, 
the diagram gives the Fi and F 2 form factors as functions of the form Fi_ 2 (g^/'^^, p/m, e). 
The form factor can be expanded as a power series in g^/m^ at hxed e, followed by the limit 
6 ^ 0 , 


Fi = F, ( 0 ) 


^ ^ + (A„ + Bo)\og^ + D„ 

^UV ^IR ^ 

, 2 — + — + {Ai + Bi) log — + Di 

dq^ leuv (^ir rn 


2 /r\\ 

+? — (0) 


+ . 


( 12 ) 
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FIG. 1. Non-abelian contribution to the one-loop vertex correction. 


and similarly for the F 2 form factor. It is conventional to label e as either euv or 
depending on whether the integral is nltraviolet or infrared divergent. Ultraviolet divergences 
are cancelled by renormalization connterterms. Infrared divergences cancel when a physically 
measnrable process is compnted. Expanding the form factor in and then taking the limit 
e —0 gives an expression that is analytic in and misses terms which are non-analytic in q^. 
The non-analytic terms are not needed for the calcnlation of the coefficients in the effective 
theory, since the effective Lagrangian is analytic in momentnm. The coefficients of the 
effective Lagrangian are determined by compnting (for example) the difference of Fi in the 
fnll theory and effective theories. The non-analytic terms in Fi cancel in the difference, and 
the analytic terms determine the nnknown parameters cp ■ ■ ■ cb 2 in the effective Lagrangian. 

Loop diagrams in HQET are fnnctions F({p},p,e) times powers of the coefficients q 
in the effective Lagrangian, where {p} are the external momenta. All on-shell loop graphs 
vanish when expanded in powers of p, followed by e —> 0. This is becanse the coefficient of 
any power of p is a dimensionally regnlated integral of the form 



There is no dimensionfnl parameter in the integrand, so the integral vanishes. The matching 
condition is then trivial: one takes Eq. (JT^), and throws away the 1/e terms to obtain the 
difference of the graph in the fnll and effective theory. All the 1/e terms in the difference are 
nltraviolet divergences (which are cancelled by renormalization connterterms), since there 
are no infrared divergences in matching conditions. To see this more explicitly, one can 
evalnate integrals snch as Eq. (^) by breaking them np into the snm of two integrals, one 




only ultraviolet divergent, and the other only infrared divergent. For example, 


f k 1 

f k 

1 m^ 

1 

1 

r 1 

1 1 

/ (27r)'^A:4 

J (27r)‘^ 

k^ {k? + m^) (fc^ -|- m^) 

IbTT^ 

■ ^UV 

^IR- 


since euv = ^ir = e. A given quantity in the effective theory is of the form 


( 14 ) 


A,n{- -!-V (18) 

\^UV ^IR/ 

There can be no hnite parts (the analogs of the {A + B) log/i/m and D terms in Eq. ([l^ ), 
since the net integral is zero. A typical matching condition is of the form 


graphs in full theory = graphs in effective theory + q. 


(16) 


where q is a coefficient in the effective Lagrangian. Using Eq. (|T^ and Eq. ([T5|) , the 
matching condition can be written as 


- \- {A + B) log — + D — —Aeff-h Cj. 

eiR m eiR 


(17) 


The 1/euv terms are cancelled by the renormalization counterterms in the full and effective 
theories, respectively. The coefficients in the effective Lagrangian have no infrared diver¬ 
gences. Thus B = — Aeff, and 




Ci — {A-\- B) log- \~ D, 


m 


(18) 


i.e. Ci is obtained from Eq. (|T^) by keeping the hnite pieces, and omitting the 1/euv and 
1/eiR terms. 

The coefficients of the /2m) Q and (D^/8m^) Q are hxed by the dispersion 

relation of QCD, C 2 = C 4 = 1. The other terms, which all contain at least 

one power of the gauge held A^, are obtained by computing the one-loop Q^QA on-shell 
scattering amplitude. The wave-function renormalization graph Fig. ^ and the vertex cor¬ 
rection Fig. 1^ can be found in many textbooks on quantum held theory |^. In dimensional 
regularization, one hnds that the wave function graph is 
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FIG. 2. One-loop wavefunction renormalization. 
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where 


FIG. 3. Abelian one-loop vertex correction. 
iS (p) = -iC 2 iQ) ^ [a m + B (p^) i)) 


A (p^) = J dx T (e/2) (4 — e) m^x — p^x(l — x) ^ 

B (p^) = ~J dx T {e/2) (2 — e) {1 — x) m^x — p‘^x{l — x) ^ 


C 2 (Q) = T“T“ = 


(19) 

( 20 ) 
( 21 ) 


is the Casimir of the quark representation. The on-shell wavefunction renormalization cor¬ 
rection is 


= -6*2 (g) 


Att 


B R) + 2 m= {A + 


= c's (g) 


o 


TT 


1 1 ^ 3 , m 

- -^-h 1 — - log — 

2e;7y eiR L p 


( 22 ) 


The on-shell vertex correction Fig. ^ can be expressed in terms of the form-factors Fi 
and F 2 , 


igT^ u (p') 


Fi{q^) 7^ + *^ 2 (g') 


2 \ qi' 


2m 


(p). 


(23) 


where q = p' — p. The form factors are 

b'P (9O = A fc2 (Q) - Ic, (ad) 
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and 


where 


F2^'’ {q^) = ^ (<^2 (Q) - ^C2 (ad)^ I . 

I (q^/m?) = f dx — --r 

^ 'Jo — q^x[l — X) 

= I'dx 

J 0 w 


K {q^= J dx 


m 


rn? — q^x{l — x) 


log 


rn? — q^x{l — x) 
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(26) 

(26) 

(27) 

(28) 


and 


C 2 {ad) = 3 

is the Casimir of the adjoint representation. Expanding to order q^/n? gives 


Pi'"' = ^ (Cj (<3) - \c-- (‘“I)) 

fP = ^ (C 2 (Q) - icj (ad)) 


+ 


1 ^ 3. m g^/ll 1 1, m' 


2e[/y ^iR 
1 

- + 


+ 1 - - log-h 




-1— log ■ 

3e,H 8 3 ^ ft 


_2 12 m 2 _ 

The hnal diagram is the non-abelian vertex correction Fig. ^ This is compnted in 
backgronnd held Feynman gange, which preserves gange invariance. The resnlting diagram 
can also be evalnated in terms of the Fi and F 2 form-factors, 
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The total on-shell form factors at one loop are given by 


F 2 


1-5Z + 


( 33 ) 
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71 m? 
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3e,R 8 3 


m 




C 2 {ad) 


^{V) 


<9) 


(34) 




TT 




as q 


TT 


The total form-factor Fi (0) is unity, since gauge invariance is preserved by the background 
held method. 

The scattering amplitude for a low-momentum heavy quark off a background vector po¬ 
tential can be computed by expanding Eq. and multiplying by ^Jm/E for the incoming 
and outgoing quarks. If p is the three-momentum of the incoming quark, p' is the three- 
momentum of the outgoing quark, and q = p^ — p, one hnds that the effective interaction 

is 


Unr (pO 
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unr (P) 


(35) 


where 









(36) 


and 


j =i (p + p') + i''<1 - if (IpI'+ ip'i V <1 


z 

2m 


16m^ 
+F 2 (g') 
1 

16m^ 


z 

8m^ 
1 


(37) 


pf - Ipt) - X (p + p') - ^ (Ipf + IPl'O (P + (Ip' 


-CT X q- - 

2m 16m'^ 


'|2 I |2\ 

p I - IpI q- 


8m^ 

i i i2 

q| CT X q 


2 I |2\ 

- IpI q 


^ I 1 2 / / . \ 

q| (p +p) 


16m^ 


8m^ 


(IpI^ - IpI^) CT X (p' + p) + ^cr ■ (p' + p) (p' X p) y 


11 




























Comparing Eqs. (^,^) with the scattering amplitude in the effective theory from the 
Lagrangian Eq. gives 


cp — Fi + F 2 (38) 

cd = F^ + 2 F 2 + SF[ (39) 

cs = Fi + 2 F 2 (40) 

cwi = -^1 + 2-^2 + 4F( + AF 2 (41) 

cw2 = \f2 + AF[ + AF!2 (42) 

Cp'p = F 2 (43) 

cm = -If2-4FI (44) 


where 


p = F,(o), f; = 


dF 


(46) 


(} 2=0 


d /m?) 

Note that the nine parameters (including C 2 and C 4 ) in the effective Lagrangian Eq. (Q) are 
determined in terms of only three independent constants, E 2 , F[, and F 2 , since Fi = 1. 


Reparametrization invariance [jl^ gives six linear relations among the coefficients. This will 
be discussed in more detail in the next section. 

The explicit expressions for the coefficients are obtained using Eqs. 
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(47) 

(48) 

(49) 
(60) 
( 61 ) 
(62) 
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FIG. 4. Quark contribution to the vacuum polarization. 
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FIG. 5. Gluon contribution to the vacuum polarization. 


The results for NRQED can be obtained by setting C 2 (ad) = 0 and C 2 (Q) = 1, and agree 
with those found in with the replacement 


5 

log/r ^ log2A - - 
o 


(53) 


The difference in finite parts is because the NRQED integrals were evaluated in Ref. 
using a momentum space cutoff, instead of using dimensional regularization. The results for 


the 1/m operators agree with known results for HQET The 1/m^ matching conditions 


at tree-level, and the /i dependence at one-loop also agree with known results Note 

that cp is independent of /i in QED. This is easy to see if one computes the renormalization 
of the magnetic moment operator in the effective theory in Coulomb gauge, in which all 
transverse photon interactions are suppressed by 1/m. 

The discussion so far has concentrated on the fermion part of the effective Lagrangian. 
There is, in addition, the pure gauge held part of the effective action. The one-loop correction 
to the gluon propagator is shown in Figs. (^) and (^. The gluon diagram is the same in QCD 
and in HQET, so the one-loop matching condition is from the quark vacuum polarization 
diagram. This gives the effective action il5|-|l7 



with 
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( 55 ) 


a 


di = l- (Q) logmV/i^ 
Stt 


a 


where 




T(Q) = 1/2 


is the index of the quark representation. The identity 


0 = / 2 B" G-"'' + 2 9 Ubc G;)„Gf,G° + G;!„ G-"'' 


has been used to eliminate Dy from Eq. 


V. REPARAMETRIZATION INVARIANCE 


The coefficients of operators in the HQET Lagrangian are constrained by reparametriza- 
tion invariance [M. The reparametrization invariant spinor held is given by 




(56) 


where 'ip^ is the conventional heavy quark held that satishes 


= 'ip^, 


(57) 


A{w,v) is the Lorentz transformation matrix 


A {w, v) = 


1 + pvp 


2 (1 + w • n) 


and 


= 


+ iD^/m 


( 68 ) 


(59) 


\v^^ + iDf^/ml 

One needs to choose a particular operator ordering for the covariant derivatives; diherent 
orderings are related to each other by held redehnitions. 
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It is simplest to consider the consequences of reparametrization invariance when 
in Eq. (^). Then there is no operator ordering ambiguity, and the held can be written 
as 


= 


1 1/2 


2\J {m + id ■ v)^ + {id^y + id ■ v + \J{m + id ■ v)^ + 


(60) 


X 


m + +id ■ V + i^^ + \J{m + id ■ v)^ + {idi_y 




where 


d^ = d^ -v^ d-v. 


(61) 


If one uses Eq. (|60|) , replaces tjjy by the spinor UNRe~^^'^, with j^unr = unr, 

u^nr^nr = 1, and v = (1,0, 0,0)), the held reduces to the spinor that satis- 

hes the Dirac equation '^u = mu, and is normalized so that uu = 1. This shows that 
reparametrization invariance will determine the coefficients of the 1/m suppressed operators 
which are hxed by relativistic invariance. 

The reparametrization invariant kinetic term is 




{m + id ■ v)^ + {id±f — m 


i’v 


(62) 


This is not the same as the terms in the Lagrangian Eq. (|^). The reparametrization invariant 
held Eq. (^) does not automatically produce a Lagrangian in the “standard” NRQCD form. 
However, one can convert Eq. 


ijv = 


to this form by making a held redehnition, 

1/2 


{m + id ■ V) + {idi_) + m 
_m + id ■ V + — {id^y + m^ 

The kinetic energy term in the primed held is 






m + id ■ V — \J — {idi_)^ + m^ 


<■ 


(63) 


(64) 


which when expanded gives Eq. (^ with C2 = C4 = 1. Thus C2 = C4 = 1 follows from 
reparametrization invariance. The transformation factor in Eq. (|5B[ ) when applied to on- 
shell spinors (instead of helds) reduces to E. This is the same as the hux factor for the 
incoming and outgoing particles which was included in Eq. (|35D . 
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To determine the constraints of reparametrization invariance on the effective Lagrangian, 
consider Eq. (^) with the gauge helds included, i.e. with d ^ D. Expanding to order 1/m^ 
gives 




V 


1 + /1 + 


2m 


i^v 


(65) 


where 


A = 1 - 


{iv ■ D) 




+ 


4m^ 


B - I- ^ ^ • D)^ 


m 


8m‘^ 




( 66 ) 

3). A different ordering 


A particular ordering has been chosen for the operators in Eq. (^ 
gives an effective Lagrangian that is related by a held redehnition. 

The most general reparametrization invariant Lagrangian is a linear combination of in¬ 
variant terms, such as + Hj)/m}j ^Gvi etc. The effective Lagrangian ob¬ 

tained in this way is not in the form Eq. (P), but it can be converted into that form by held 
redehnitions that preserve 'fi'ipv = i’v One hnds by a straightforward (but not very enlight¬ 
ening computation) that the ehective Lagrangian is a linear combination of the invariant 
linear combinations 


iv ■ D O2 + O4 + Op -\- Op) -\- Os + Oy\/i, 

20p 40/) -|- 4 O 5 -|- Owi + Ow2 + ‘^Opip — Omi 

Owi + Ow2, ( 67 ) 

20/) -|- Owi + Ow2 — Omi 
0/11, Oa 2, Oyi3, Oa4, Obi, Ob2, 

up to terms of order 1/m"^. Here O/?, etc. are the operator coefficients of c/?, etc. in Eq. (||). 
The above linear combinations imply the constraints 

C2 = 1, 

C 4 = 1, 
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cs = 2 cf — 1 , 


( 68 ) 


CW2 — Cwi — 1 , 

Cp'p Cp 1, 

2cm = cf — cd- 

which are satished by Eqs. (^),(^^- 

VI. CONCLUSIONS 

The HQET/NRQCD Lagrangian has been computed to one loop and order 1/m^, and 
has been shown to be reparametrization invariant to order a/m?. The original form of 
the NRQCD propagator Eq. (|) cannot be used to compute the effective Lagrangian by 
matching to QCD. Instead, one must treat the propagator as an inhnite series, and resum 
the series after doing the loop integral. As a result, the matching computations for NRQCD 
and HQET are the same. It is straightforward to obtain the effective Lagrangian (in the 
one-quark sector) to higher orders in 1/m by expanding the form factors Ti 2 and the spinors 
in the computation of Eqs. (^) , (^6|) to higher orders. No further Feynman graphs need to 
be evaluated. 
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